Dynamic Structure in a Molten Binary Alloy by ab-initio Molecular Dynamics: 
Crossover from Hydrodynamics to the Microscopic Regime. 
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The dynamic structure factor of the ^Lio.eiNao.sg Uquid alloy at T = 590 K has been calculated 
by ab initio molecular dynamics simulations using 2000 particles. For small wavevectors, 0.15 < 
q/ < 1.6, we find clear side peaks in the partial dynamic structure factors. Whereas for q < 0.25 
the peak frequencies correspond to the hydrodynamic sound dispersion of the binary alloy, 
for greater g-values we obtain two modes with phase velocities above and below the hydrodynamic 
sound. A smooth transition between hydrodynamic sound and the two collective modes is shown 
to take place in the range 0.25 < q/ < 0.35. The mass ratio in this system, mNa/m-Li ~ 3, is 
the smallest one so far for which the fast mode is observed. We also predict that inelastic X-ray 
scattering experiments would be able to detect the slow mode, and explain why the inelastic neutron 
scattering measurements [P.R. Gartrell-Mills et al, Physica B 154, 1 (1988)] do not show any of 
these modes. 

PACS numbers: 61.20.Ja, 61.20.Lc, 61.25.Mv 
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Interest in the collective excitations in liquid binary 
mixtures was generated by the pioneering molecular dy- 
namics (MD) study on the dynamics of the liquid Na-K 
alloy 111 . It was further stimulated by the MD results for 
liquid Li4Pb, where a new, high-frequency mode, sup- 
ported by the light Li atoms only (the so-called "fast 
sound") was identified [Q. Kinetic theory calculations 
1^ for two-component fluids with large atomic mass dif- 
ference, confirmed the appearance of a collective mode 
at high frequencies, lu, and wave numbers q > qh, where 
qu denotes the upper limit of the hydrodynamic regime. 
Moreover, this mode propagated with a phase velocity 
close to that of the pure light component, which was 
clearly greater than the hydrodynamic sound velocity, Ch, 
of the two-component system. Subsequently, calculations 
based on the revised Enskog theory for binary mixtures 
of hard spheres with a large mass difference Q predicted 
two propagating collective modes in addition to hydro- 
dynamic sound. One of these was identified with the 
fast sound while the other mode, which had propagat- 
ing phase velocities below Ch, was consequently termed 
"slow sound". Theoretical |Q, ^, 0, computer simula- 
tion @, g I, |l0| and experimental f|, 0, |l|, |l||l|, |l| 
studies have investigated the existence and properties of 
these collective modes. Experimental evidence for the 
fast sound has been obtained for systems such as Il2-Ar, 
Il2-Xe and He-Xe mixtures by light scattering [|l5|, and 
He-Ne, He-Ar and Li4Pb by inelastic neutron scattering 
(INS) 1^, p. Note that all these systems except Li4Pb 
are gas mixtures, whereas the latter is pseudoionic mix- 
ture, and that the mass ratio varies from around 33 for 
He-Xe to around 5 for He-Ne. 

There has been discussion concerning the way in which 
the two collective modes connect with the hydrodynamic 
one when qh is approached from above. Some theoretical 
and MD results for He-Ar and He-Ne mixtures have sug- 



gested that either the fast and slow sound merge into the 
hydrodynamic sound [0, ||] , or that the fast sound disap- 
pears at qh and the slow sound merges into the hydrody- 
namic one, with a predicted value ||, @, || for qh ~ 0.07 
A^^. A similar value for q^, has recently been obtained 
for liquid Li4Pb by both INS and by MD calculations 
1^, |l^. Recent MD calculations |p^, w hich followed the 
INS measurements of Bafile et al |l4[|, for the dynamic 
structure of a Heo.77Neo.23 gas mixture at two densities, 
showed a clear crossover from hydrodynamics to fast and 
slow modes at 0.2 < q/ < 0.5 and gave qu ~ 0.2, 
which is substantially greater than previous estimates. 

The present communication addresses the previous 
questions while providing further insight into the micro- 
scopic dynamics of the Lio.6iNao.39 liquid alloy. This 
is a metallic system, with a rather small mass ratio « 
3, to which much theoretical and experimental attention 
has been devoted because of its strong phase separating 
tendencies. It has a consolute point at T « 577 K and 
concentration CLi ~ 0.64, close to the "zero alloy" com- 
position (cLi ~ 0.61) for which the total static structure 
factor, S{q), reduces to the concentration partial struc- 
ture factor Scc{q), i-e. S{q) = Scc{q)/cuCN&, because 
of the negative scattering length of the isotope ^Li. The 
S{q) has been measured at several temperatures and con- 
centrations jl^, and INS measurements ||l^ have been 
performed for several temperatures at the zero alloy com- 
position. The measured total dynamic structure factor, 
S{q,bj), in the explored {q^uj) region decreased monoton- 
ically as a function of to, and side peaks, which would 
have been indicative of collective modes, were absent. 

The present study has used the orbital-free ah ini- 
tio molecular dynamics (OF-AIMD) method, where the 
forces acting on the nuclei are computed from electronic 
structure calculations, based on density functional theory 
(DFT), which are performed as the MD trajectory is gen- 
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erated. A simple liquid metallic alloy, AcBi_c, is treated 
as a disordered array of bare ions interacting with the va- 
lence electrons through electron-ion potentials. The total 
potential energy of the system can be written, within 
the Born-Oppenheimer approximation, as the sum of 
the direct ion-ion coulombic interaction energy and the 
ground state energy of the electronic system, Eg[pg{f))\. 
According to DFT, the ground state electronic density, 
Pg{r) minimizes an energy functional which is given as 
the sum of the kinetic energy of independent electrons, 
Ts[p], the classical Hartree electrostatic energy, EH[p]^ 
the exchange-correlation energy, E^c[p\, for which we 
have adopted the local density approximation and finally, 
the electron-ion interaction energy, i?ext[p], for which 
we have used local ionic pseudopotentials constructed 
within DFT In the OF-AIMD approach Q ||^ 

an explicit but approximate functional of the density is 
used for Ts[p]. Proposed functionals consist of the von 
Weizsacker term, Tw[p{r)\ = |/ c^^l Vp(r)p/p(r), plus 
further terms chosen in order to reproduce correctly some 
exactly known limits. Here, we have used a simplified ver- 
sion of the average density models, developed by Garci'a- 
Gonzalez et al pTf in which Tg = T\y + Tp, where 

Tp-yJ drpir^/'-'^Hrr (1) 
fc(f) = {2k°pf j dsk{s)wp{2k°p\r-s\) 

k{r) — (Stt^)^/'^ the Fermi wavevector for 

mean electron density poi and Wf3{x) is a weight function 
chosen so that both the linear response theory and the 
Thomas- Fermi limits are correctly recovered. Further de- 
tails of this functional are given in reference and we 
merely note that in the present simulations we have used 
/3 = 0.51. 

We have performed OF-AIMD simulations for the 
Lio.6iNao.39 liquid alloy at temperature T = 590 K and 
number density p= 0.03218 A^"^. The cubic simulation 
box contained 2000 particles. Given the ionic positions at 
time t, the electronic density is expanded in plane waves, 
and the energy functional is minimized with respect to 
the plane wave coefficients yielding the ground state elec- 
tronic density and energy. The Hellmann-Feynman theo- 
rem is used to obtain the forces on the ions which, along 
with the Verlet leapfrog algorithm, are used to update the 
ionic positions and velocities. The timestep was 0.0025 
ps, the equilibration lasted for 10 ps and the calculation 
of properties was made averaging over another 60 ps. The 
cutoff energy for the plane wave expansion was 8.15 Ryd, 
giving a basis of about 85000 planewaves. These choices 
allow a minimum g- value of 0.158 which will permit 
investigation of length and time scales covering a range 
which includes the hydrodynamic regime. 

The total dynamic structure factor 5^(9,0;), which is 
directly related to the intensity obtained in either an in- 



elastic neutron (INS) or X-ray scattering (IXS) experi- 
ment, is a weighted average of the partials Sij{q,uj), de- 
fined as 

g».(g.^) = 2^(^/^^.)i/2 / dte^-HpMt)-p-{qM 

(2) 

where Ni is the number of i-type particles, the (....) 
stands for the ensemble average, the asterisk denotes 
complex conjugation, pi{q,t) = X]/(i)=i exp[ig • Ri{i}{t)] 
represents the density fluctuations of the i-th component 
with wavevector q , and (<) is the position of the 
i-type particle I at time t. The 5t((7, observed in an 
IXS experiment is given by 

S^-'\q,.) . ± (c.c,r/^^^5,(<z,.) (3) 

where fi{q) are the atomic scattering factors and {P{q)) 
=^^^i Ci fiiq)- Correspondingly, in an INS experiment 
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where the S^{q,u;) are the self-dynamic structure factors, 
(bi) is the coherent scattering length, 47r(6^) is the total 
scattering cross section and (&^) = '^ii^i) is the 

average cross section per atom. 

The results for Suuiq,^) and 5'NaNa('3', w) are shown 
in fig. |l| for three representative g- values. At the smallest 
q- value allowed by the simulations, q= 0.158 A~^, both 
SuLi{q,(^) and 5NaNa(9,<^) show clear Brillouin peaks at 
a common frequency, ujb, which is the typical behaviour 
in the hydrodynamic region and indicates a propagating 
sound mode with a hydrodynamic velocity Ch = wb / q 
« 3000 m/s. 5'LiLi('?,w) exhibits Brillouin peaks for 
the other two g-values of fig. |^ and continues to show 
side peaks or shoulders up to q w 1.5 A~^, whereas 
S^aNaiq, has only a shoulder for q — 0.48 A"-^ and this 
feature disappears at q « 0.9 A~^. In contrast, the par- 
tial longitudinal currents, Cij{q,uj)= uP' S'ij(q,a;), show 
clear peaks at any g-value. At the lowest g-value, the 
frequencies of the maxima in the Cij (<?, w) coincide with 
those of the Brillouin peaks, but at larger q they are 
somewhat higher because of the contribution from the 
damping terms 0. Moreover, as q increases C\i{q,iS) 
diminishes pointing to a progressive decoupling of the 
motions of two species. 

The dispersion curves WLiLil?) and WNaNal?) of the in- 
elastic peaks in 5'LiLi('3', w) and 5'NaNa('7, i^) respectively. 



3 




10 20 ' 10 20 30 ■ 20 40 



mips'') 

FIG. 1: OF-AIMD results, at three g-values, for the hq- 
uid Lio.6iNao.39 at T=590 K. Upper frame: Partial dynamic 
structure factors (logarithmic scale). SLiLi{q,(j->) (full lines) 
and 5'NaNa('Z,w) (dashed lines). Lower frames: Partial longi- 
tudinal currents: Chihi{q,io) (full lines), CNaNa(iJ, t^;) (dashed 
lines) and CLiNa('?,<^) (dot-dashed lines). 



are shown in fig. ||. Up to q ~ 0.3 A""'^, ujuuiq) and 
'^NaNa('z) show the Same linear behaviour of hydrody- 
namic collective excitations propagating with speed Ch 
« 3000 m/s. Above this g- value, the dispersion curve 
splits into two branches signaUing the onset of dynamic 
decouphng of the species, and giving rise to fast and slow 
sound modes in the alloy. The fast mode involves the Li 
particles only, with a phase velocity Cfast =^uu{q) /<1 ~ 
3800 m/s, whereas the slow mode phase velocity is Csiow 
= ijJNaNa(9) / 1 ~ 2100 m/s, at q—O.&b A^^. We have 
also performed OF-AIMD simulations for pure Li at the 
same temperature and total number density as the alloy, 
and the calculated sound velocity, 3900 m/s, is very close 
to the velocity of the fast mode. Following usual prac- 
tice 0, we have also obtained dispersion curves derived 
from the maxima in Cij{q,uj), and these are also shown 
in fig. ^. The frequencies for both Li and Na are larger 
than LUuu{q) and a;NaNa('?), and the earlier onset of the 
branching suggests that the departure from hydrodynam- 
ics begins already at g « 0.15 A~^. Moreover, the higher 
slope obtained for the Li leads to a larger estimation for 
the phase velocity of the fast mode (cfast ~ 4600 m/s). 

Now we consider the possible experimental observa- 
tion of these collective modes. As already mentioned, 
the INS measurements |lj] gave a total structure factor 
with no side peaks as a function of ui. We have calculated 
, to) using eqn. (^) and data from our simulations. 
The results illustrated in fig. ^ do not show side peaks 
in agreement with the experiment, and despite the ap- 
pearance of collective modes in the simulation. For an 
understanding of this anomaly we have also plotted in 
fig. ||the partials Sij{q,oj). We see that 5'LiNa(9,w) also 
exhibits a peak at a frequency between that of S'LiLi(<Z, ^) 
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FIG. 2: OF-AIMD results for the dispersion curves of the 
collective modes in 5'LiLi(?,'^) and S'NaNa(g, tt') (full circles 
and squares respectively). The positions of the maxima of 
CLiLi(g, ti.') and CNaNa(g,t<j) are also shown (open circles and 
squares respectively). 



and S'NaNaC'Ziw). In the weighted average of eqn. (||), the 
negative scattering length of the ''Li isotope ((&Li)=-0.22 
X 10^^^ cm) gives rise to a negative contribution from the 
partial S'LiNa(9,'^) which balances the contributions from 
the peaks in S'LiLi(9, and S'NaNaC?, w). In contrast, we 
predict that side peaks would be observed in an IXS mea- 
surement because the atomic form factors involved in the 
average in eqn. (^ always take positive values. Moreover, 
the peaks would be those due to the slow sound because 
the prefactor of S^aX^J^q^u:) is about twenty times larger 
than that of S'LiLi(<Z, "j^) and twice that of 5LiNa('?, 
with the overall result that the the total Sj~ {q,uj) is 
mainly controlled by S'NaNa('3', i-^). 

In conclusion, our simulations of the Lio.6iNao.39 liq- 
uid alloy show well defined collective modes analogous 
to those already found in other binary fluids. However, 
this system has the smallest mass ratio considered so far 
and is the first metallic system for which "fast sound" 
appears. The results for the INS total dynamic structure 
factor explain the reason for the failure to detect collec- 
tive excitations in the INS experiment by Gartrell-Mills 
et al [r^ . However, we predict that these excitations, 
particularly the slow sound mode, would be observed in 
an IXS experiment. 

In order to specify where the transition from hydro- 
dynamics to kinetic regime takes place, we have found 
that it is important to specify which magnitude is used. 
In particular, partial currents depart from hydrodynamic 
behavior for smaller q- values than the corresponding par- 
tial structure factors, and this must be taken into account 
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FIG. 3: OF-AIMD results for the total dynamic structure 
factors in liquid Lio.6iNao.39 at T=590 K. Open circles: 

Si^'^^\q,u)). Full circles: si^'^^\q,uj). For comparison, the 
partials Shiu{q,to) (full line), SNaNa(?,t<^) (dashed line) and 
>S'LiNa(9,w) (dotted line) are also included. 



in discussions of the onset of "fast sound" . 

Finally, we have shown that is now possible to per- 
form ab-initio molecular dynamics simulations to study 
the dynamic properties of fluids. The OF-AIMD method, 
which employs the electron density as the basic variable, 
gives an approximate treatment of the electron kinetic 
energy, but within this limitation allows the simulation 
of large samples for long runs. These are basic require- 
ments for performing reliably the ensemble averages that 
define most dynamic properties in liquids, as well as for 
studying the ordering tendencies in multicomponent sys- 
tems. 
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